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Transverse momentum integrated multiplicities in the central region of pp collisions at LHC ener- 
gies satisfy Koba-Nielsen-Olesen scaling. We attempt to relate this finding to multiplicity distribu- 
tions of soft gluons. KNO scaling emerges if the effective theory describing color charge fluctuations 
at a scale on the order of the saturation momentum is approximately Gaussian. From an evolution 
equation for quantum corrections which includes both saturation as well as fluctuations we find that 
evolution with the QCD /3-function satisfies KNO scaling while fixed-coupling evolution does not. 
Thus, non-linear saturation effects and running-coupling evolution are both required in order to 
reproduce geometric scaling of the DIS cross section and KNO scaling of virtual dipoles in a hadron 
wave function. 
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I. INTRODUCTION 

The color fields of hadrons boosted to the light cone are thought to grow very strong, parametrically of order 
A + ~ 1/g where g is the coupling [l| . The fields of nuclei are enhanced further by the high density of valence charges 
per unit transverse area, which is proportional to the thickness A 1 / 3 of a nucleus 0. 

In collisions of such strong color fields a large number of soft gluons is released. Due to the genuinely non- 
perturbative dynamics of the strong color fields a semi- hard "saturation scale" Q s emerges; it corresponds to the 
transverse momentum where the phase space density of produced gluons is of order l/a s . The mean multiplicity per 
unit rapidity in high-energy collisions is then n = (dN/dy) ~ l/a s . Below we argue that a semi-classical effective 
theory of valence color charge fluctuations predicts that the variance of the multiplicity distribution is of order 
k^ 1 ~ 0{oP s ) so that the perturbative expansion of n/k begins at order l/a s ^> 1. We show that in the strong field 
limit then a Gaussian effective theory leads to Koba-Nielsen-Olesen (KNO) scaling Q. This relates the emergence of 
KNO scaling in p^-integrated multiplicity distributions from high-energy collisions to properties of soft gluons around 
the saturation scale. 
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FIG. 1: KNO scaling of charged particle multiplicity distributions in non-single diffractive pp / pp collisions at various energies 
as measured by the UA5 Q, ALICE [g and CMS [g] collaborations, respectively. Note that we restrict to the bulk of the 
distributions up to 3.5 times the mean multiplicity. 



The KNO scaling conjecture refers to the fact that the particle multiplicity distribution in high-energy hadronic 
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collisions is universal (i.e., energy independent) if expressed in terms of the fractional multiplicity z = n/n. This is 
satisfied to a good approximation in the central (pseudo-) rapidity region at center of mass energies of 900 GeV and 
above 0, Q as shown in fig. [TJ On the other hand, UA5 data [1] taken at y/s = 200 GeV appears to show a slightly 
distorted multiplicity distribution. This is in line with the observation that at lower energies higher-order factorial 
moments G q of the distribution are energy dependent and significantly different from the reduced moments C q [7| : 

^ (n (n-l)---(n-q+l)) _ 

In fact, since the difference of G q and C q is subleading in the density of valence charges one may interpret this 
finding to indicate that the high density approximation is less accurate for y/s = 200 GeV pp collisions. Approximate 
KNO scaling has been predicted to persist also for min-bias p + Pb collisions (at LHC energies) in spite of additional 
Glauber fluctuations of the number of participants and binary collisions [1] . A more detailed discussion of multiplicity 
distributions at TeV energies is given in refs. and references therein. 

Transverse momentum integrated multiplicities in inelastic hadronic collisions are not governed by an external 
hard scale, unlike say multiplicity distributions in e + e~ annihilation or in jets (loj . Hence, the explanation for the 
experimental observation should relate to properties of the distribution of produced gluons around the saturation 
scale Q s . 

II. KNO SCALING FROM A GAUSSIAN ACTION IN THE CLASSICAL LIMIT 

We shall first discuss the multiplicity distribution of small- x gluons obtained from a Gaussian effective action for 
the color charge fluctuations of the valence charge densities p Q , 

Z = J Vpe' SMv[p] , (2) 

S MV [p] = J d xx J_Jx . (3) 

In the strong field limit a semi-classical approximation is appropriate and the soft gluon field (in covariant gauge) can 
be obtained in the Weizsacker- Williams approximation as 

A+(Z-, X± ) = -g-Lp a ( Z -, X± ) =g J d 2 Z ± J -p°(«- Z X ). (4) 

Parametrically, the mean multiplicity obtained from the action is then 

n Q s b±_ , (5) 

a s 

where Sj_ denotes a transverse area and Q s ~ g 2 p. The prefactor in (J5J can be determined numerically [ill [l2[ but 
is not required for our present considerations. 

One can similarly calculate the probability to produce q particles by considering fully connected diagrams with 
q valence sources p in the amplitude and q sources p* in the conjugate amplitude (for both projectile and target, 
respectively). These can be expressed as [1J| 1 



d q N 



dyi---dy q/ iiimi 

where the reduced moments 



)con„. Cq (dyi) {dy q ) 



(6) 



This expression is valid with logarithmic accuracy and was derived under the assumption that all transverse momentum 
integrals over pt.\ • • 'PT.q are effectively cut off in the infrared at a scale ~ Q s due to non-linear effects. 



1 The rapidities y\ ■ ■ ■ y q of the q particles should be similar. Here we assume that all particles are in the same rapidity bin. 
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The fluctuation parameter k in eq. ([7]) is of order 

fc~(i\£-l)#Sx 



(8) 



Once again, the precise numerical prefactor (in the classical approximation) has been determined by a numerical 
computation to all orders in the valence charge density p [15(. 

The multiplicity distribution is therefore a negative binomial distribution (NBD) [l3j], 



P(n) 



T(k + n) 



n n k k 



T(k)T(n + 1) (n + k) n + k 



(9) 



Indeed, multiplicity distributions observed in high-energy pp collisions (in the central region) can be described quite 
well by a NBD, see for example refs. [lij. The parameter fc _1 determines the variance of the distribution 2 and can 
be obtained from the (inclusive) double-gluon multiplicity: 



d 2 N 
dyidyi 



dN 
dyi 



dN\ 
dy2/ 



(10) 



From this expression it is straightforward to see that the perturbative expansion of fc _1 starts at 0(a®) since the 
connected diagrams on the lhs of eq. (|10[) involve the same number of sources and vertices as the disconnected diagrams 
on the rhs of that equation (also see appendix). This observation is important since in general the NBD ^ exhibits 
KNO scaling only when n/k 3> 1, and if k is not strongly energy dependent. A numerical analysis of the multiplicity 
distribution at 2360 GeV, for example, achieves a good fit to the data for n/k ~ 6 — 7 |8|, which we confirm below. 
Such values for n/k have also been found for peripheral collisions of heavy ions from ab initio solutions of the classical 
Yang-Mills equations [lq ]; furthermore those solutions predict that n/k < 1 for central collisions of A ~ 200 nuclei. 

To illustrate how deviations from KNO scaling arise it is instructive to consider a "deformed" theory with an 
additional contribution to the quadratic action. We shall add a quartic operator [l7| . 



Sq[p] 



d 2 v ± 



- fp a 0i -vi)p"K ,vj_) 



2/2 2 (t^ 



dvn 



p a (v 1 ,v ± )p a (v 1 ,v ± )p b (w 2 ,v ± )p b (u 2 ,v_l) 



h'4 



fii) 



We assume that the contribution from the quartic operator is a small perturbation since K4 ~ gigA 1 ' 3 ) 3 while 
A 2 ~ 9 (d A 1 / 3 )- In the classical approximation the mean multiplicity is unaffected by the correction as it involves 
only two-point functions 3 . On the other hand, k^ 1 as defined in (|10[) now becomes 



AL 2 



2tt 



1 - 3/3 (iV c 2 + 1) 



with j3 



Therefore, in the classical approximation 



n N r , , 



6ir 3 Qi 



1)) 



«4 



dz p A {z ) 



(12) 



(13) 



This result illustrates that n/k decreases as the contribution of the ~ p A operator increases. We repeat that the 
derivation assumed that the correction is small so that (fT3"|) does not apply for large values of (3N 2 . 

Ref. [l7| estimated by entirely different considerations that for protons (3 ~ 0.01 at x = 10~ 2 . That would 
correspond to a smaller value of n/k by a factor of 1.43 than for the Gaussian theory. Assuming that RG flow with 
energy approaches a Gaussian action [1S| . n/k should increase by about this factor. NBD fits to the data shown in 
fig. [2] confirm that n/k indeed increases with energy, which might indicate flow towards a Gaussian action; however, 
the observed increase from y/s — 200 GeV to 7 TeV is much stronger: a factor of about 3. This apparent discrepancy 
could be resolved at least partially by running of the coupling in eqs. (|5I13[) with Q s but this requires more careful 
analysis 4 . 



2 More precisely, the width is given by n Vfc -1 + ~ fi/y/k; the latter approximation applies in the limit n/k 3> 1, see below. 

3 The two-point functions (pp) in the theories J3} and need to be matched. Thus, the "bare" parameters p? in (|3]l and p? in 
are different as the latter absorbs some self-energy corrections. We refer to ref. |17( for details. 

4 In fact, the running of a s at the effective scale Q s is taken into account if the mean multiplicity is computed with energy evolved 
unintegrated gluon distributions like e.g. in refs. [a. Il2ll . 
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FIG. 2: KNO scaling plots of charged particle multiplicity distributions at \rj\ < 0.5 in NSD collisions at various energies and 
NBD fits; z = N c h/{N c h) and ^{z) = (N c h) P(N c h). Note that the mean multiplicity quoted for the fits has been rescaled by 
1.5 to include neutral particles; also, that here k is integrated over the transverse plane of the collision. 



III. QUANTUM EVOLUTION AND THE DISTRIBUTION OF DIPOLES IN THE HADRONIC WAVE 

FUNCTION 

In the previous section we considered the multiplicity distribution of "produced gluons" in a collision of classical YM 
fields sourced by classical color charges p moving on the light-cone. At high energies though (i.e., when a s log x^ 1 ~ 1) 
the classical fields are modified by quantum fluctuations [l[ . Resummation of boost-invariant quantum fluctuations 
leads to an energy dependent saturation scale, for example, as required in order to reproduce the growth of the multi- 
plicity n with energy. In particular, the energy dependence of the mean saturation scale, averaged over all "evolution 
ladders" (distribution of quantum emissions) , can be obtained by solving the running-coupling BK equation [l9| . 

Instead, in this section we shall solve a stochastic evolution equation which accounts both for saturation (non-linear) 
effects as well as for the fluctuations of the rapidities and transverse momenta of the virtual gluons in the wave function 
of a hadron before the collision. We do this in order to determine the multiplicity distribution (rather than just the 
mean number) of dipoles in a hadronic wave function boosted to rapidity Y. 

We shall do so by solving via Monte-Carlo techniques the following evolution equation for P[n(x),Y], which is the 
probability for the dipole size distribution n(x) to occur: 



f z [n(x) - S xz ] P[n(x) - S XZ ,Y] - / f z [n(x)} P[n(x), Y] . (14) 



dP[n{x),Y] 
dY = 

Note that in this section x = log 1/r 2 denotes the logarithmic dipole size (conjugate to its transverse momentum) 
rather than to a light-cone momentum fraction. This equation has been studied before in ref. 20] for fixed a s and in 
ref. [2ll ] for running a s (r 2 ). Those papers also provide references to related earlier work. 

The first term in (|14p is a gain term due to dipole splitting while the second term corresponds to loss due to 
"recombination" . 

fz[n{x) = j-— (15) 



T z [n{x)] = 1 - exp / n(x) log (1 - t(z\x)) (16) 



is the splitting rate and 



is the dipole scattering amplitude for a dipole projectile of size z to scatter off the target with the dipole distribution 
n(x). Note that T z [n(x)] is non-linear in the dipole density as it involves also the pair (and higher) densities. Finally, 

r 2 

r(x\y) = a(x)a(y) exp(-|a; - y\) = a(r<) a(r 2 ) (17) 

r > 
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is the elementary dipole-dipole scattering amplitude at LO in perturbative QCD. For more details we refer to ref. [21|. 
Here, we recall only that it was found there that evolution with a running' coupling suppresses fluctuations in the tails 
of the travelling waves and so restores approximate geometric scaling |22j |. 

We have determined the multiplicity distribution of dipoles with size ~ 1/Q S (Y), 



Ni(Y) 



dx rii(x, Y) 



(i = 1---10 5 ) 



(18) 



r 2 <l/(Q*(Y)) 



by evolving a given initial configuration n(x, Y — 0) 10 5 times. Despite starting with a fixed initial condition, evolution 
introduces fluctuations in the rapidities where splittings occur, and in the sizes of the emerging dipoles. 
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FIG. 3: Multiplicity distribution of virtual dipoles with size r 2 < 1/(Q 2 S (Y)) in the wave function. Left: evolution with trivial 
/3-function, a s =const. Right: QCD /3-function. 

In fig. |3] we show that fixed-coupling evolution does not obey KNO scaling of the distribution of virtual quanta while 
running-coupling evolution does. The shape of the distribution however looks different than the measured distribution 
of produced particles from fig. [T] This could be due to the fact that our evolution model does not treat diffusion in 
impact parameter space. Hence, P(N) shown in fig. [3] should be interpreted as the multiplicity distribution at the 
center of the hadron. 
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Appendix A: The moment C2 = k 1 with the quartic action 



We can obtain the fluctuation parameter k by calculating the inclusive double gluon multiplicity and expressing it 
in terms of the single inclusive or mean multiplicity. The connected two particle production cross section for gluons 
with rapidity y\ and y 2 has the form: 

d 2 N \ /dN\/dN\_ / d 2 N 

dN* 



N 2 (p,q) \ dyidy2 / \dy 1 /\dy 2 ) \ dy x dy 2 / co „ n . ^ 
is the mean multiplicity and the brackets denote an average over events. N 2 {p, q) is given by: 
12 r — d 2 k t L tl {p 1 k 1 )L^p,k 2 )L„(q,k 3 )L^q,k 4 ) 

X 



9 f 

N 2 (p,q) = J^^fgaa'fg'bb'fgcc'fg'dd' J 



^ r ga aJ go DJ g CCJ gaa j 11 ( 27r ) 2fc a (p _ fcl )a(p _ ^ _ ^ _ ^ 

(pr(k2)pl\k i )p<i(k 1 )p d 1 (k 3 ))(pf(p- k 2 )pf (q- fc 4 W(p- ^{q- k 3 )) . 
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L M denotes the Lipatov vertex, for which: 



4fc 2 

L^p,k)L»{p,k) = r (p-k) 



(A2) 



For the four-point function in the target and projectile fields we use [TtJ 

(pf(P - k 2 )pf{q - h) P C 2(p - h)pi'(q - k 3 )) = 



(27T) 



dz jl 2 {z ) 



6 a ' b 'd c ' d 'S(p + q-k 2 ~ k 4 )S{p + q-k 1 ~k 3 ) 



-6 ac 'S t " 1 S(k 1 - k 2 )S(k 3 - k 4 ) + 6 a ' d '6 b ' c 'd{ P ~q~k 2 

l 2 



-(2tt)< 



dz pL (z ) 



(xa'b'ftc'd' +S a'c' 6 b'd' 



.ga'd'gb'c' 



h)S(p — q — ki 
S(k 1 + k 3 - k 2 



The first two lines on the rhs of the above equation originate from the quadratic part of the action while the third line 
is due to the quartic operator. The product of the Gaussian parts of the two four-point functions gives nine terms, 
one of which (~ 5 ac 5 bd S a c S b d ) corresponds to a disconnected contribution. It exactly cancels the second term in 
eq. (|ATj) . 
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FIG. 4: One of eight connected diagrams for two-gluon production with the quadratic MV action. 

Four of the other eight terms (~ S ac 5 bd or - 6 a ' c 'd b ' d ') give identical leading contributions to double gluon pro- 
duction. They correspond to a "rainbow" diagram like the one shown in Fig. 01 In the "rainbow" diagram, on one 
side (target or projectile), the p's corresponding to the same gluon momentum are contracted with each other. The 
remaining four "non-rainbow" diagrams are suppressed relative to the terms we keep at large p and q [l3| . Hence, 
the leading Gaussian contribution is: 



.9 

8tt 7 



dz fi 2 {z ) 



S ± N 2 (N 2 - 1) 



4 4 

p q 



(A3) 



The same reasoning applies also for the additional quartic contribution and only "rainbow" diagrams are considered, 
like the one in Fig. [5] There are two of them (one for the projectile and one for the target) to first order in 1 , and 
their contribution is: 



9 f f t f ff\ * 

2(2tv) g aa cc / (2n 

i—l 



d 2 h L„(p, W(P> k 2 )L„{q, k 3 )L»{q 7 ki) 



) 2 fc 2 (p - fci) 2 (p - k 2 Y{q - k 3 f{q - fc 4 ) 2 



2(2tt) 8 

7T 2 K4 



dz fi 2 {z ) 
(s a ' b 's c ' 



dz fi (z ) 



ra'c' rb'd' 



S^d^Sih - k 2 )S{k 3 - k 4 ) x 



ca'd' xb'c' 



S(ki + k 3 - k 2 - fc 4 ) 



The color factor evaluates to 



FIG. 5: Connected diagram for two-gluon production from the quartic operator in the action ]1_ 



f g aa>f g -wf g cc>f g - dd >8 ac 5 bd (s a ' b 'S c ' d ' + S a ' c '6 b ' d ' + 6 a ' d '5 b ' c ') = 2N 2 (N 2 - 1) + N 2 C {N 2 C - if 
Using eq. (|A2|) we get: 



(27r) 8 7T 2 K 4 



dz fx (z ) 



dz fl (z ) 



[2N 2 (N 2 -l) + N 2 (N 2 -l) 2 ] x 



S± 



d 2 fci 



d 2 k. 



p 2 q 2 J k 2 (p—ki) 2 J k 2 (q — k 3 ) 2 
The integral over the ladder momentum is again cut off at the saturation scale Q s : 

d 2 k\ 2ir p 

klip-krf ~ p og qJ • 

Then, the quartic contribution to connected two gluon production becomes 



<J 



12 



n 2 



dz fi 2 (z ) 



1 2 



dz fi A (z ) 



Si 



[2N 2 (N 2 C - 1) + iV c 2 (iV c 2 - I) 2 ] log -f log ■ 

P Q s 



47T 8 K 4 

The last step is to express the fully connected diagrams in terms of the single inclusive cross section 

P 



dz fx (z ) 



4tt 4 



Summing eq. (IA3[) and eq. (|A4[) and using eq. (|A5|) we get 



I dN 
\dy 



d 2 N 
dy\dy 2 



2n 



4(iV c 2 + l) [Jdz-^(z-)]'- 

Q 2 {N 2 - 1)S X KiS±(N* - 1) [fdz-piz-)]' 



/dN\/dN 
\ dyi I \ dy 2 



(A4) 



(A5) 



The fluctuation parameter k 1 is now identified with the expression in the square brackets. We rewrite it in terms of 



Cp 9 8 
6ir 3 Q 2 K4 



and use 



to arrive at the final expression 



Q 



2 _ 9 4 Cf 
2tt 



dz fi^(z ) 



dz fx (z ) 



iV 2 -l o2 „ 1 



Qis x ~ = 1 - mN 2 c + 1) 

2n k 



(A6) 



(A7) 



(A8) 
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Appendix B: The moment C3 with the quartic action 



In this section we are going to calculate the connected diagrams for three-gluon production to obtain the correction 
to the reduced moment C 3 at order 1/k 4 ~ l/[g(gA 1 ^ 3 ) 3 ] : assuming as before that gA 1 ^ 3 > 1. At the end of this 
section we also outline corrections suppressed by higher powers of gA 1 ^ 3 . 

We are looking for the contribution of the connected diagrams to the following expression [23| : 



d 3 N 



dyidy 2 dy 3 



8(2tt 

6 

n 



ig 3 f 

q~ fgaa 1 fg'bb' fg"cc' fgff fg'ee' fg"dd' ^ 



d 2 h L a {p, fci)L Q (p, k 2 )L (q, k 3 )LP(q, k 4 )L, ( (l, k 5 )L~<(l, k 6 ) 



LI (2n) 2 k 2 (p_ kl f(p- fca)2(g- k 3 ) 2 (q- fc 4 ) 2 a - W- fee) 2 
(p^(p - k 2 )pT{q - k 4 )pl d (l - k 6 )pU P - - fc 3 )p?(« - *B)> X 



(fc 2 )P2 e {k A )pf (k 6 )p a 2 (fc x )p§ (fc 3 )p§ (fe 5 )) 



(Bl) 



As before, the p correlators of the target and the projectile consist of two parts, one from the quadratic operator in 
the action and another from the additional p 4 operator: 

{p p p p p p ) = {p p p p p p /Gaussian + \P P P P P P /Correction ■ 

The product of the two Gaussian contributions from the target and the projectile, to leading order in the gluon 
momenta, gives rise to 16 "rainbow" diagrams. The result has been obtained previously jl 31 ] and reads (expressed in 
terms of the mean multiplicity) : 



d 3 N 



8tt 2 



dN\/dN\/dN\ 



QiS\{N 2 - l) 2 \d Vl /\ dy 2 /\dy 3 / 



x dyidy 2 dy 3 / Conn Gaussian 
The correction, to first order in kJ 1 is 

~ 2 \P P P P P P /Gaussian \P P P P P P /Correct 



(B2) 



(B3) 



Again, we are considering only rainbow diagrams, so for the Gaussian six-point function in the above expression, from 
all possible contractions, we keep only the term 



dz p. (z ) 



8 a ' s '5 h ' e '8 c ' d '8{k 1 ~ k 2 )5(k 3 - k 4 )S(k 5 - k 6 ) 



To calculate the correction to the six-point function to first order in k 4 we factorize it into a product of two- and 
four-point functions. There are fifteen possible factorizations of that kind. Three of them are disconnected diagrams 
and the remaining twelve give identical contributions. We consider, for example, the following combination: 

(pl f (p - k 2 ) P T{q - h)pf(l - h)pl(p - kxjpHq ~ ka)pl(l - fc 5 )> 

= (pUp - h)pl(q - k 3 ))(pl f (p - k 2 )p?(q - k 4 )pl d (l - h)p\{l - fc 5 )) . 

The two point function is 



(pKp-k 1 )p b 1 (q-k 3 )) = (2ir) i 



dz p, (z ) 



5 ab 5(p - 



and for the correction to the four-point function we use the last line from eq. (|A3[) . 
The color factor is 

fgaa'fg'bb'fg"ee'fgff'fg'ee'fg"dd'S ab (S^ + ST 8* + ^ '«J de ) ? 5*' e ' S C ' d ' 



2N 3 (N 2 



N?(N 2 C 
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Putting everything together into eq. (jBll) and multiplying by two [because of (|B3|) ] and by twelve (which is the number 
of possible diagrams) we get: 



d 3 N 



dy\dy 2 dy 3 



is. 



3 / Conn. Corrcctic 



K4 





4 




J dz p 2 (z ) 




J dz p 4 (z ) 



[2iV c 3 (7V c 2 _ i) + Ar3(iV c 2 - l) 2 ] x 



d 2 fc t L a (p, fci)L a (p, k 2 )Lp(q, k 3 )L^{q, fc 4 )L 7 (Z, fc 5 )£ 7 (i, fa) 

L = l (27T)2fc2 (p - fc x )2(p _ fc 2 )2( g _ fc 3 )2( g _ fc 4 )2 (Z _ ^(J _ fc 6 )2 X 

5(ki - k 2 )S(k 3 - k i )S(k 5 - k 6 )S(p + q - k x - k 3 )S(p + q - k 2 - fc 4 - k 6 + k 5 ) 



n 



3g 18 S± 
16K47T 13 p 2 q 2 l 2 





4 




J dz p 2 (z ) 




J dz p 4 (z ) 



[2N*{N*-l) + N?(N*-lf] x 



d 2 k x 



fc?(p + g-fci) a (p-*i) 4 J k\{l-k 2 ) 



d 2 k 2 



Again, we regularize the ladder integrals at the saturation scale, 

d 2 k 2tt 1 



/ 



k 2 (p + q-k) 2 (p-k) 4 p 2 q 2 Q 2 ' 
Finally, using expression (IA5|) for the mean multiplicity the p contribution to three-gluon production becomes 

d 3 N \ _ 487r(iV 2 + l) [/ dz~p 4 (z^)] 2 / dN\ / dN\ / dN 

. Correction ~ ~ K 4 Q 2 s Sj_(N* - l) 2 [J dz - p {z - j] 2 \dy~J \dy~ 2 / \dy~ 3 



Summing (JB2]) and (|B4|) . 



(B4) 



d 3 7V 



dyidy 2 dy 3 / G 



8tt 2 



4847V 2 + 1) [Jrfz-^(z-)] 1 

QfS 2 (iV 2 - l) 2 « 4 Q 2 ^i(iV 2 - l) 2 [/ds-/i 2 («-)] : 



/ dN\ / dN\ / dN\ 
\dyi/\ dy 2 J \ dy 3 / 



From the above equation the third reduced moment is: 



C 3 = 



8tt 2 



or 



487r(iV 2 + l) [Jdz-p 4 (z-)Y 

QjSl(N? - l) 2 Ki Q 2 S 2 ± (N 2 - l) 2 [fte-pb-ff 



M_ll!g4 5i ^ = 1 _ 9 ^ (iVc2 + 1) 



(B5) 



where we have used expressions (|A6|) and (|A7[) . 

For a NBD we have that C 3 = 2/k 2 but if we compare (|B5|) to the square of eq. (|A8[) . which is 



4tt 2 



(B6) 



we see that the coefficients of the corrections at order O{0) differ. That means that the p 4 operator in the action 
provides a correction to the negative binomial distribution. 

In fact, such deviation from a NBD is more obvious if even higher order operators are added to the action. Dropping 
the longitudinal dependence of the operators for simplicity, such an action would have the form 



S ~ / d 2 vj. 



§ab p a p b d abc p a p b p c S ab S cd + perm. a b c d S ab d cde + perm. a b c d e 



2p 2 



K3 



K4 



-p p p p 



«5 



-p p p p p 



(S ab S cd 6 ef + perm.) + <d abc d def + perm.) „ , r d e f 

'-p a p b p c p d p e p t + 



K6 



The additional terms are suppressed by powers of gA 1 ! 3 (24|: 

p 2 ~ g{gA^ 3 ) , n 3 ~ g{gA^ 3 ) 2 , Ki ~ g{gA^ 3 ) 3 , n 5 ~ g{gA^ 3 ) 4 , n e ~ j^ 1 / 3 ) 5 
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The cubic operator gives a correction to the six-point function, i.e. to C3 at order but does not correct the 

four-point function, i.e. Ci = \jh (it only renormalizes fj, 2 ). The same applies to the p 6 operator: C3 will contain a 
term ~ but 1/k does not. Hence, beyond a quadratic action the relation C3 = 2/k 2 is not exact. 
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